Disorder-Induced Order in Two-Component Bose-Einstein Condensates 
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We propose and analyze a general mechanism of disorder-induced order in two-component Bose- 
Einstein condensates, analogous to corresponding effects established for XY spin models. We show 
that a random Raman coupling induces a relative phase of 7r/2 between the two BECs and that the 
effect is robust. We demonstrate it in ID, 2D and 3D at T = and present evidence that it persists 
at small T > 0. Applications to phase control in ultracold spinor condensates are discussed. 
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Degenerate quantum gases offer unprecedented control 
tools opening fascinating possibilities, e.g. investiga- 
tions of quantum disordered systems ■ Current ex- 
perimental 0, 0, 0| and theoretical [171, [3, Hi works 
are mainly devoted to studies of the interplay between 
disorder and nonlinearity in Bose-Einstein condensates 
(BEC) in a quest for traces of Anderson localization. In 
the regime of strong correlations, evidence for the Bosc 
glass phase has been reported [1], and even more exotic 
quantum phases have been proposed [2, [lo| . 

Weak disorder can have strong effects also in classical 
systems. For instance, a general mechanism of random- 
field-induced order (RFIO) has been proposed recently 
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^ ^ It is responsible for ordering in graphene quan- 
tum Hall ferromagnets [l3l, "^He-A in aerogel and amor- 
phous ferromagnets [l^ , as well as for inducing superflu- 
idity in hardcore bosonic systems [111 ■ This effect is best 
understood in classical ferromagnetic XY models in the 
presence of uniaxial random magnetic fields. For the 2D 
square lattice, the Hamiltonian reads 



H 



i'-ji=i 



(1) 



where the spins are unit 2D vectors in the XY plane: 
CTi — {cos 9i, sin 9i) at site i € Z^. When h = the sys- 
tem does not magnetize as a consequence of the Mermin- 
Wagner-Hohenberg (MWH) theorem [ill. In contrast, 
a weak uniaxial random field breaks the continuous 
U{1) symmetry. Then, the MWH theorem does not ap- 
ply and the system spontaneously magnetizes with a non- 
zero component of the magnetization perpendicular to 
the random field. This has been proven at zero tempera- 
ture and strong arguments have been given that the effect 

Hamiltonian ([1]) can 
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persists at small temperatures 
be realized with ultracold atoms in optical lattices, 
the effect is rather weak 11 1. 



but 



In this Letter, we propose an analogue of the RFIO ef- 
fect using two BECs trapped in harmonic potentials and 
coupled via a real- valued random Raman field. We show 
that the mean field Hamiltonian of the two-component 
BEC is analogous to the XY spin Hamiltonian ([T]), with 
the Raman coupling playing the same role as the mag- 
netic field in Eq. ([T]), and the relative phase between 
the BECs corresponding to the spin angle 6i. Then, the 
RFIO effect shows up in the form of a relative phase be- 
tween the BECs fixed at a value of ±tt/2. The finite-size 
two-component BEC system is continuous and formally 
equivalent to the discrete spin system ([T]) on an infinite 
lattice. We find that even in low dimensions, the RFIO 
effect is much more pronounced and robust in coupled 
trapped BECs than it is in uniform lattice spin models. 
Note that trapped (finite size) BECs at sufficiently small 
T show true long range order also in ID and 2D as phase 
fiuctuations take place on a scale larger than the size of 
the systems [l^l- We demonstrate the effect in ID, 2D 
and 3D at T = and present strong evidence that it 
persists for small T > 0. 

Interestingly, the RFIO effect is quite general. For 
instance, consider the two-spin lattice Hamiltonian: 



H = - 



E 



(2) 



where 17^ arc independent real-valued random couplings 
with (identical) symmetric distributions. In this system, 
it can be proven rigorously that there is no first order 
phase transition with the order parameter ct,; • in di- 
mensions d < 4 [ill. More precisely, in every infinite- 
dimensional Gibbs state (phase), the disorder average of 
the thermal mean (crj • Tj) takes the same value. By sym- 
metry, this value has to be zero, implying that the av- 
erage cosine of the angle between ai and Ti is zero. At 
T = 0, these results also apply [ill and arc consistent, 
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by analogy, with the relative phase 7r/2 of two randomly 
coupled BECs, discussed below. 

We consider a trapped two-component Bose gas with 
repulsive interactions and assume that the two compo- 
nents consist of the same atomic species in two different 
internal states, coupled via a position-dependent (ran- 
dom, quasi-random, or just oscillating) real-valued Ra- 
man field fl{r) of mean zero (/ fidr = 0). The typi- 
cal amplitude and spatial variation scale of 51 (r) are de- 
noted by and Ar . At sufficiently small T, the trapped 
gases form BECs which can be represented by the clas- 
sical fields V'i.2(r) in the mean-field approximation. The 
energy functional of the system then reads 



E = J dr[ {h^/2m)\Vi:i\ 



V^(r)|^i|2 + (51/2)1^-1 



+ (?iV2m)|VV'2|' + ^(r)|^2p + (.92/2)1^21^ 

+giM^\H' + m{r)/2) (V'i>2 + V'2>i)] , (3) 

where V(r) is the confining potential, and g,; = ^■nh^ai/m 
and gi2 = 47rft^ai2/m are the intra- and inter-state cou- 
pling constants, with ai and ai2 the scattering lengths 
and m the atomic mass. The last term in Eq. ([3)) repre- 
sents the Raman coupling which can change the internal 
state of the atoms. 

The ground state of the coupled two-component BEC 
system is obtained by minimizing E' as a function of the 
fields V'l and il}2 under the constraint of a fixed total 
number of atoms N = J dr(|-0ip -|- |'02p)- This leads to 
a set of two coupled Gross-Pitaevskii equations (GPE): 

+ mm^J, (4) 

with /i the chemical potential and i = 2(f) for i = f (2). 

At equilibrium, for fl^ = and gi,g2 > gi2, the 
BECs are misciblc [l^. Their phases 9i are uniform, 
arbitrary and independent. Now, a weak Raman cou- 
pling (SjriRl ^ /i) does not noticeably affect the densi- 
ties. However, arbitrarily small fl{r) breaks the continu- 
ous U{1) symmetry with respect to the relative pha se of 
the BECs and, following the results of Refs. [Ill, [H, Ell , 
the relative phase can be expected to be fixed. To 
make this clearer, we neglect the changes of the densities 
when the weak Raman coupling is turned on and ana- 
lyze the phases. For simplicity we suppose gi = 52 and 



p(r) = pi(r) = p2(r)- The substitution 4'i 



in the energy functional ^ leads to E = Eo + AE where 
Eq is the energy for SIr = and 



AE 



drp{r) 



r n2 

—(wef + hn{r) cose 

Am 



drp(r)|l(Ve)^ 



(5) 



where Q = 9i + 62 and 9 = 9i — 92- Minimizing AE im- 
plies = const; hence the second line in Eq. (O vanishes 



and the only remaining dynamical variable in the model 
is the relative phase 9 between the BECs. Note that if 
Pi 7^ P2 the variables Q and 9 are coupled and one cannot 
consider them independent (the pi ^ p2 case is analyzed 
in the sequel). Equation ([5]) is equivalent to the classical 
field description of the spin model H]) in the continuous 
limit, where the relative phase 9{v) represents the spin 
angle and the Raman coupling ri(r) plays the role of the 
magnetic field. Thus, we expect RFIO [lH to show up 
in the form cos 9 ~ for weak random r2(r). 

Let us examine Eq. ([5]) in more detail. It represents a 
competition between the kinetic term which is minimal 
for uniform 9 and the potential term which is minimal 
when the sign of cos 61 is opposite to that of 51 (r). For 
Ml-^x. 3> h^/2mX^, the potential term dominates and 9 
will vary strongly on a length scale of the order of Ar. In 
contrast, if ftf^R ^ h? /2mX^ the kinetic term is impor- 
tant and forbids large modulations of 9 on scales of Ar. 
The Euler-Lagrange equation of the functional ([5|) is 



2m 

V [p(r)V6l] + — p(r)r2(r) sin6i 



0. 



(6) 



For the homogeneous case (p ~ const) and for slowly 
varying densities (neglecting the term Vp), assuming 
small variations of the relative phase, 9{r) = 6*0 + S9{r) 
with \d9\ <C TT, the solution of Eq. © reads 

d9{k) ~{2m/h){n(k)/\k\^)sm9o (7) 
in Fourier space. Inserting Eq. ([7|) into Eq. ([5]), we find 



AE^ -mp / dk (|r2(k)|V|k|^)sin2 6io. 



(8) 



The energy is thus minimal for 9q — ±7r/2, i.e. cos 9q — 0. 
This indicates RFIO in the two-component BEC system 
owing to the breaking of the continuous U{1) symmetry 
of the coupled GPEs. For a random Raman coupling, 
even if the resulting fluctuations of 9 are not small, the 
average phase is locked at 6*0 — ±7r/2. Note that if 9{r) 
is a solution of Eq. so is —9{r). This follows from 
the fact that for any solution (■0i,V'2) of the GPEs 
(i/'i , f/'l) is also a solution with the same chemical poten- 
tial. The sign of 6*0 thus depends on the realization of 
the BECs and is determined by spontaneous breaking of 
the 9^—9 symmetry. 

Let us turn to numerics starting with gi ~ g2- For 
homogeneous (p = const) gases, we solve Eq. (|6]). Fig- 
ure [1] shows an example for a ID two-component BEC, 
where n{x) is a quasi-random function chosen as a sum 
of two sine functions with incommensurate spatial peri- 
ods. The dynamical system ^ is not integrable. It turns 
out that the solution we are interested in corresponds to 
a hyperbolic periodic orbit surrounded by a considerable 
chaotic sea. Figure[T]eonfirms that 9{x) oscillates around 
9q ~ ±7r/2. The oscillations of 9{x) are weak and fol- 
low the prediction ([7]), which in ID, after inverse Fourier 
transform, corresponds to the double integral of Q{x). 
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Figure 1: RFIO effect in a ID two-component BEC 
trapped in a box of length L and in a quasi-random Ra- 
man field. Panel (a); Raman coupling function f2(a;) = 
-100(;i/2mL2)[sin(2;/AR + 0.31) + sin(a;/(2.44AR) 1.88)] 
with Ar = 0.00939L. Panel (b): Relative phase e[x) = 
Qiix) — 02 (s) obtained by solving Eq. ^ numerically (solid 
black line) and comparison with Eq. ((7)l (dashed red line — 
nearly identical to the solid black line). 



For trapped gases and for g\ ^ g2 we directly solve the 
coupled GPEs Figure [2] shows the results for a ID 
two-component BEC in the Thomas-Fermi regime con- 
fined in a harmonic trap with a random fl(x). A typical 
realization is shown in Fig. [^tt. For each realization of 
0(a;), the resulting relative phase 6 can change signifi- 
cantly but only on a scale much larger than Ar because 
/ir^R <^ h'^/2mX'^, as shown in Fig. [2b- However, averag- 
ing over many realizations of the random Raman coupling 
and keeping only those with / 6{x)dx > (resp. < 0), 
we obtain {0{x)) w 7r/2 (resp. — 7r/2), with the standard 
deviation about O.Stt as shown in Fig. 

The dynamical stability of the solutions of the 
GPEs (|3]) found in the ID trapped geometry can be 
tested by means of the Bogoliubov-de Gennes (BdG) the- 
ory which allows also to estimate the quantum depletion 
of the BECs The BdG analysis shows that the solu- 
tions of the GPEs ([4]) are indeed stable and that the BdG 
spectrum is not significantly affected by the Raman cou- 
pling. It implies that turning on the Raman field does 
not change the thermodynamical properties of the sys- 
tem, and the RFIO effect should persist for sufficiently 
low T > 0. Note that the GPEs (g]) possess also a solu- 
tion with both components real. However, this solution 
is dynamically unstable. In fact, there is a BdG mode 
associated with an imaginary eigenvalue and the corre- 
sponding BECs phases (under a small perturbation) will 
evolve exponentially in time. In addition, the BdG anal- 
ysis shows that the quantum depiction is about 1% and 
can therefore be neglected. 

Calculations in 2D and 3D — whose detailed results 
will be published soon — show essentially the same 



Figure 2: RFIO effect in very elongated (effectively ID) 
trapped BECs. The data corresponds to *^Rb atoms in two 
different internal states in an anisotropic harmonic trap with 
frequencies uj^ = 2tt x 10 Hz and uj± — 2Tr x 1.8 kHz. The 
total number of atoms is TV = 10^ and the scattering lengths 
are ai = 5.77 nm, 02 — 6.13 nm and ai2 = 5.53 nm. Panel 
(a): Single realization of the random Raman coupling hO./ ^ 
for Ar = 10"^Ltf and ft^R ~ 3 x 10"^/i. Panel (b): Relative 
phase 9 corresponding to ^{x) shown in panel (a). Panel (c): 
6 averaged over many realizations of Q,{x) (solid line) and the 
averaged value ± standard deviation (dashed lines). In panel 
(c) the solutions with J 6dx > only are collected (the other 
class of solutions with 6 —9 is not included). 



Disorder- Induced Ordering effect in all dimensions. For 
example, Fig.[3]shows the result for two coupled 3D BECs 
in a spherically symmetric harmonic trap. Here, the Ra- 
man coupling is a sum of quasi-random functions similar 
to that used for Fig. [T] in each spatial direction and with 
Ml^ ~ 10~^/i. The density modulations arc found to 
be negligible. However, even for this low value of the 
Raman coupling. Fig. [3] shows that the relative phase is 
fixed around 9q = 7r/2 with small fluctuations. Other 
calculations confirm that the sign of 6q is random but 
with 1 6*0 1 = 7r/2 for all realizations of VI{t) and that the 
weaker the Raman coupling, the smaller the modulations 
of 6{v) around 6^. This shows once again the enormous 
robustness of RFIO in two-component BECs. 

In summary, we have shown that RFIO occurs in a 
system of two BECs coupled via a real-valued random 
Raman field. It has been demonstrated in ID, 2D and 
3D for homogeneous or trapped BECs. The signature of 
RFIO is a fixed relative phase between the BECs around 
9{) = ±7r/2. For quasi-random Raman coupling, the fluc- 
tuations can be very small (O.OStt for the parameters used 
in Fig. [1]). For completely random Raman coupling the 
fluctuations can be larger (about 0.37r for the parameters 
used in Fig. Interestingly, the two-component BEC 
system is continuous and RFIO is stronger and more ro- 
bust than in lattice spin Hamiltonians of realistic sizes 
llj. RFIO can thus be obtained in current experiments 
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Figure 3: RFIO effect in a 3D two-component BEC trapped 
in a spherically symmetric harmonic trap with frequency to = 
27r X 30Hz. The total number of atoms is A'' = IC , the scatter- 
ing lengths are as in Fig.[5]and we use a quasi-random Raman 
coupling Sl(a;, J/, 2;) oc E„e{a;,H,z} [sin('i/AR) +sin(u/(1.71AR))] 
with Ar = 4.68/im/27r and ftfiR ~ 5 x 10~''/i. Shown is the 
relative phase Q in the plane z = 0/^m in units of tt. 



with two-component BECs j2ll . l22l | and observed using 
matterwave interferometry techniques 22 1. 

Apart from its fundamental importance, RFIO can 
have applications for engineering and manipulations 
of quantum states by providing a simple and robust 
method to control phases in ultracold gases. We find 
particularly interesting applications of phase control in 
spinor BECs and, more generally, in ultracold spinor 
gases [l|. For example, in a ferromagnetic spinor 
BEC with = 1 as in ^^Rb, the wavefunction is ^ oc 
(e-*-^ cos2(6i/2), V2sin(6»/2) cos(6'/2), e+"^ cos2(6'/2)), 
the components correspond to mp = 1,0,-1 
and the direction of magnetization is n — 
(sin Q cos (/>, sin d sin 0, cos ff) . Applying two real- 
valued random Raman couplings between mp — Q and 
mp = ±1, fixes </> = or tt, i.e. the magnetization 
will be in the XZ plane. By applying two random 
real-valued Raman couplings between nip = and 
rap = \ and between rap = —\ and rap = 1, we force 
the magnetization to be along ±Z. Similar effects 
occur in antiferromagnetic spinor BECs with F = 1, as 
^^Na. Using Raman transitions with arbitrary phases, 
employing more couplings, and higher spins F offers a 
variety of control tools in ultracold spinor gases. 
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